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We have studied the magnetic and thermodynamic properties as well as the NMR spectra of 
the quasi two-dimensional Heisenberg antiferromagnet RbFe(MoC>4)2. The observed temperature 
dependence of the order parameter, the critical indices and the overall magnetic H — T phase 
diagram are all in a good agreement with the theoretical predictions for a 2D XY model. The 
temperature dependence of the specific heat at low temperature demonstrates a crossover from a 
T 2 law characteristic of a two-dimensional antiferromagnet to a three-dimensional T 3 law. 

PACS numbers: 75.50.Ee; 76.60-k. 



I. INTRODUCTION 



The problem of finding the ground state of a two- 
dimensional triangular lattice antiferromagnet (TLAM) 
is of particular interest due to the possibility of find- 
ing solutions with unconventional magnetic order, which 
are influenced significantly both by frustration and by 
zero-point fluctuations. As the magnetic exchange in- 
teractions between the ions located on a regular two- 
dimensional triangular lattice are frustrated, a usual 
Neel-type ground state with anti-parallel alignment of 
the nearest neighbour moments is hindered by the geom- 
etry. In the molecular field approximation, the minimum 
of energy is reached for a planar spin configuration with 
the magnetisation of the three sublattices arranged at 
120° to one other. Such a spin configuration has been 
found in numerous compounds with a stacked triangu- 
lar lattice^, regardless of the dimensionality of the mag- 
netic interactions, which can have a pronounced ID or 
2D character. A very recent report suggests that instead 
of a long range structure, a spin-liquid state with a rela- 
tively short correlation length may be stabilised in some 
compounds 2 . 

In an applied magnetic field, the mean-field ground 
state of the 2D triangular antiferromagnet remains highly 
degenerate, as the overall energy is defined only by the 
total spin of the three sublattices, while there is an in- 
finitely large number of different states with the same 
total spin^i. Indeed, any possible value of the total 
spin can be represented as an umbrella-type structure 
of equally tilted magnetic moments with identical trans- 
verse components, or, alternatively, as a variety of both 



planar and nonplanar structures with variable angles be- 
tween the magnetic moments and the direction of the 
field. 

This degeneracy is removed by thermal and quantum 
fluctuations, which select a symmetric planar structure 
for a purely isotropic (Heisenberg) case, for an easy-axis 
type of anisotropy and also for an easy-plane type of 
anisotropy if the field is applied in the basal plane 5,6,7,8 . 
If the field is applied perpendicular to the easy-plane, 
the mean-field approach yields an umbrella-type struc- 
ture, which is non-degenerate with respect to the mutual 
orientation of magnetic moments. For the planar struc- 
tures, fluctuation analysis reveals a characteristic feature 
of the magnetisation curve M(H). Namely, a collincar 
magnetic structure with nearly one third of the satura- 
tion magnetisation (M = M s /3) that is stabilised by the 
fluctuations in the vicinity of the field H s /3, where H s is 
the saturation field. The stabilisation of this structure is 
seen as a plateau on the M(H) curve, while in the molec- 
ular field approximation the magnetisation curve has a 
continuous derivative. This structure is referred to as 
"two spins up, one down" or UUD. 

Apart from this, TLAMs are distinguished from con- 
ventional collinear ferro and antiferromagnets by the fact 
that for some phases, called chiral phases, there are still 
two ways to arrange the rest of the structure even when 
the direction of the spin on a particular apex of the tri- 
angle is fixed. Fig. ^ shows two magnetic structures, 
a and b, which share the alignment of the spin in the 
position A and possess the same exchange energy, but 
differ by the direction of rotation of the spins when trans- 
lated to the nearest apex. Chiral degeneracy influences 
considerably the nature of the phase transition to an or- 
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FIG. 1: Two different magnetic structures on a hexagonai 
iattice possessing the same exchange energy. The direction of 
magnetic moment at point A is identical for both structures. 



dered stated. Within a magnetic H — T phase diagram 
for a single TLAM, both chiral and nonchiral structures 
are possible, therefore the second order phase transitions 
with different critical indices should be expected for the 
different values of an applied magnetic field. 

An especially interesting point to consider is a tran- 
sition from a paramagnetic phase to a UUD structure 
occurring in a constant field on lowering the tempera- 
ture. This is a transition into a collinear phase with an 
uncompensated magnetic moment. Such structures with 
uncompensated magnetic moment are characteristic of 
ferrimagnets, therefore for a description of this collinear 
phase on a triangular lattice a ferrimagnetic order pa- 
rameter is used^. In ordinary ferrimagnets, however, the 
magnetic ions are located on non-equivalent crystallo- 
graphic positions, which leads to the appearance of non- 
equivalent magnetic sublattices and an uncompensated 
moment. In these conventional ferrimagnets an ordering 
transition in an applied magnetic field is absent, as an in- 
duced moment exists at all temperatures. In the case of 
a 2D TLAM, however, all the magnetic ions are in equiv- 
alent positions, while the ordering into a UUD structure 
is accompanied by a tripling of the unit cell area, which 
causes a transition of a specific type«k 

The majority of TLAMs have quasi one-dimensional 
character. These systems have been known and studied 
for some time nowi 2D TLAMs demonstrating the prop- 
erties described above are less common and have only 
recently become the subject of intense investigations. 
RbFe(Mo04)2 offers a rare opportunity to study a model 
system, which resembles closely an ideal 2D TLAM. The 
crystal structure of RbFe(Mo04)2 described by the space 
group D^ d consists of layers of magnetic Fe 3+ ions sepa- 
rated by (M0O4) 2- groups and Rb + ions. The C 3 axis 
is perpendicular to the layers. The magnetic Fe 3+ ions 
form an equilateral triangular lattice within the layers, 
whereas along the C 3 axis they are interspaced by the 
Rb + ions. Such a layered structure (depicted in Refer- 
encesll(il and lTTt) ensures the magnetic two-dimensionality 
of RbFe(Mo0 4 )2i2^. 

Relatively weak exchange interactions in 



RbFe(MoC>4)2 allow one to reach the saturation 
field experimentally, while the large spin (S = 5/2) 
of the magnetic Fe 3+ ions justify a quasi classical de- 
scription of the magnetic system. Previous experiments 
on powder and single crystal samples of RbFe(MoC>4)2 
established the presence of long range magnetic order 
below Tjv = 3.8 K and also demonstrated the stabil- 
isation of the collinear structure around M s /3***»»**«*. 
Specific heat measurements****: gave estimates of the 
magnetic entropy change during the ordering process. 
The dynamics of low-frequency excitations probed by 
ESR measurements**** gave an estimate for the ratio of 
the intraplane to interplane exchange interactions as 25 
and were also used to evaluate the magnetic anisotropy. 
By measuring the field and temperature dependence 
of the magnetisation 12 and the specific heal*****, the 
magnetic H — T phase diagram was mapped out, the 
overall shape of which agreed well with the theoretical 
predictions for a 2D TLAMiM&a. 

Neutron diffraction measurements performed in zero 
field***! have found an ordered structure with a tripled 
in-plane period, which corresponds to a 120-degree con- 
figuration within the layers and an incommensurate or- 
dering along the C 3 axis. NMR spectra of the 87 Rb ions 
located between the magnetic Fe 3+ ions from neighbour- 
ing planes were reported in ref. IT^ . These spectra have 
allowed identification of various low-temperature phases: 
a structure with incommensurate modulation along the 
axis in low fields, a tilted commensurate structure 
and a UUD structure in higher fields. 

This paper presents the results of a detailed investi- 
gation of the magnetic phase diagram and critical be- 
haviour in the different phases of RbFe(MoC>4)2 obtained 
by thermodynamic measurements. It also reports on the 
temperature dependence of the order parameter obtained 
by the 87 Rb NMR spectroscopy. The magnetic proper- 
ties, phase diagram and critical indices demonstrate good 
quantitative agreement with the theoretical predictions 
for a classical 2D XY-model. 



II. EXPERIMENTAL DETAILS 

Single-crystal samples of RbFe(MoC>4)2 were synthe- 
sised by means of spontaneous crystallisation as de- 
scribed in Ref. [1(1 The samples were thin (0.5 mm), 
nearly hexagonal plates with lengths along the edges up 
to 6 mm. The C 3 -axis was found to be perpendicular 
to the surface of the plates for all samples. The results 
reported in this paper were obtained using samples pro- 
duced by a significantly slower cool ing: of the melt com- 
pared to the samples used in Ref. Il2lll4l (1 K/hour vs 
3 K/hour). This has resulted in samples with increased 
thickness and linear dimensions. 

The magnetic moment was measured in a field of up 
to 12 T using an Oxford Instruments vibrating sample 
magnetometer. The heat-capacity measurements were 
performed by the standard heat-pulse method in the tern- 
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FIG. 2: (Color online) Magnetisation curve versus field of a 
RbFe(Mo04)2 single crystal at T = 1.55 K (top panel) and 
the derivative, 441- , taken at different temperatures (bottom 
panel). The field is applied perpendicular to the C 3 -axis. 



perature range 0.4 to 40 K in a field of up to 9 T using 
a Quantum Design PPMS calorimeter equipped with a 
3 He option. NMR spectra of the 87 Rb ions (I = 3/2, 7 = 
13.9312 MHz/T) were recorded on the home-made spec- 
trometer covering the range 35 to 114 MHz. The spin- 
echo signal was observed at constant frequency on sweep- 
ing the magnetic field in the range 2.5-9 T. Radio fre- 
quency pulses of 5 fj,s and 10 /xs were separated by the 
time interval td = 50 /is. 



III. EXPERIMENTAL RESULTS 
A. Magnetisation and specific heat measurements 

Field and temperature dependence of the magnetic mo- 
ment of RbFe(Mo04)2 are shown in Figs. [5]and|21 re- 
spectively. The M(H) and especially the dM/dH curves 
show the M s /3 plateau in the field range (H c i,H C 2) as 
well as an additional phase transition at H c3 accompa- 
nied by a weak hysteresis. The obtained data are very 
similar to the M (H) curves reported previously^, apart 
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FIG. 3: Temperature dependence of the magnetic susceptibil- 
ity of a single crystal of RbFe(Mo04)2, H _L C 3 measured in 
different magnetic fields. 



from slightly reduced hysteresis effects around the H c ^ 
field and also a small (about 0.1K) increase in the order- 
ing temperature compared to the value of Tjv obtained 
on the batch of the samples studied beforehandii2ii!ii£ 
For this new sample batch, the relative amplitude of the 
susceptibility peak at Tjy has increased by a factor of 
two for small applied fields, which usually indicates an 
improved sample quality. The anomalies, singularities 
and other special points, which are seen in the curves 
in Figs. |3 and [3] and which correspond to various phase 
transitions are used to locate the positions of the phase 
boundaries, which are later summarised in the magnetic 
phase diagram (Fig. ITT1) . 

Fig. 0] shows the temperature dependence of the spe- 
cific heat measured in RbFe(Mo04)2 for various values 
of a magnetic field applied in the basal plane of the crys- 
tal. The C(T) curves demonstrate sharp anomalies corre- 
sponding to the phase transition from the high tempera- 
ture paramagnetic phase to the low temperature magnet- 
ically ordered phase in every field from to 9 T. For the 
C(T) curves measured in the field interval 2.5 to 5.5 T, 
additional peaks are also observed when the external field 
coincides with the value of H c \. 

Fig. |S1 shows the field dependence of the specific heat 
measured at different temperatures (H _L C 3 ). For tem- 
peratures below 4 K, sharp peaks associated with the 
start of the plateau at H c \ are clearly seen on the C (H) 
curves, while the phase transitions at H C 2 or H c % could 
not be detected here. For temperatures above 4.5 K, a 
wide and round maximum in C(H) develops around a 
field of 5 T. On further heating, this maximum reduces 
in amplitude and becomes undetectable above 8 K. The 
observed peak positions were plotted on the phase dia- 
gram (Fig. Ill|) together with the data obtained from the 
magnetisation measurements. 

As can be seen from Fig. 0J the application of a weak 
magnetic field reduces marginally the ordering tempera- 
ture. However, when the field is strong enough to sta- 



4 




H H=0 T 



^ H=1T 



uH=3T 



ft, U„H=5T 



uH=7T 



uH=9T 



T(K) 



FIG. 4: Temperature dependence of the specific heat mea- 
sured in RbFe(Mo04)2, H _L C 3 measured in different mag- 
netic fields. The curves are offset along the vertical axis. 
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FIG. 5: Field dependence of the specific heat of RbFe(Mo04)2 
measured for H _L C 3 for various temperatures. 



bilise a collinear phase, the ordering temperature devel- 
ops a non-monotonic dependence with applied field. It 
rises with field for H > H c \ reaching a maximum of 4.5 K 
at H = 5 T. This maximum value exceeds significantly 
the zero-field ordering temperature of 3.9 K. 

The top panel of Fig. H3 shows the temperature de- 
pendence of the specific heat measured in RbFe(Mo04)2 
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FIG. 6: Temperature dependence of the specific heat 
(top panel) and the entropy (bottom panel) measured in 
RbFe(Mo04)2 in zero field. The dashed and solid lines cor- 
respond to the lattice and the magnetic contributions to the 
entropy. 
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FIG. 7: Temperature dependence of the specific heat of 
RbFe(Mo04)2 in the low-temperature region. A crossover 
from a square to a cubic temperature dependence is clearly 
visible on the log-log plot (the inset). 
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for a wide temperature interval in zero magnetic field. 
For temperatures in the region 10 — 40 K, the specific 
heat follows a T 2 -type of the temperature dependence. 
As this dependence is observed both below and above 
the Weiss temperature, 9cw — 25 K, it is reasonable to 
suggest that it corresponds mostly to lattice vibrations. 
Extrapolation of the lattice contribution down to zero 
temperature shows that for T < 5 K, the magnetic con- 
tribution to the specific heat is dominant and that the 
lattice contribution can be neglected. The temperature 
dependence of the entropy obtained by integrating the 
C(T)/Tvs T curves is shown in the bottom panel of Fig. 
|H] By subtracting the lattice contribution (the dashed 
line) one can obtain the temperature dependence of the 
magnetic entropy, shown in Fig. H^as a solid line. At tem- 
peratures exceeding several times the ordering tempera- 
ture Tjv the magnetic entropy tends to a value agreeing 
well with Mn(2S + 1) = TUn6. 

Fig. [7] shows temperature dependence of the zero-field 
specific heat of RbFe(Mo04)2 in the low-temperature re- 
gion. For T < 1.7 K the temperature dependence of the 
specific heat is described well by the regular cubic law 
(shown by the dash-dotted line in Fig. 01, while in the 
temperature interval 1.7 to 3.4 K it follows a quadratic 
dependence (shown by the dashed line). The inset in 
Fig. |7| reproduces the same data on a double logarithmic 
scale. 



B. NMR 87 Rb spectra 

The NMR 87 Rb spectra were recorded in the frequency 
range 35-115 MHz for various temperatures. Typical 
temperature evolution of the spectra is illustrated in 
Fig. where the top panel shows the field dependence 
of the spin-echo signal at different temperatures and the 
bottom panel shows the temperature dependence of the 
resonance fields. In the paramagnetic phase, the NMR 
spectrum consists of three lines, with the central line cor- 
responding to the (—1/2 <-> 1/2) transition and the two 
satelites to the (-3/2 <-> -1/2) and (1/2 «-> 3/2) tran- 
sitions. The frequency difference for these transitions is 
caused by the quadrupole splitting, which is found to 
be temperature independent over the entire temperature 
range under investigation. The description and inter- 
pretation of the NMR spectra in RbFe(MoC>4)2 and also 
their relation with the magnetic structure are given in our 
earlier paper—. As was shown in that paper, the value 
of the NMR frequency in RbFe(MoC>4)2 is determined 
mostly by the external magnetic field, by the quadrupole 
splitting and also by the dipolar field created by the mag- 
netic Fe 3+ ions on 87 Rb nuclei, while the contribution 
from the hyperfine field is negligibly small. 

In the paramagnetic phase, on lowering the temper- 
ature, all three NMR lines shift slightly towards higher 
fields and also broaden while approaching the ordering 
temperature T/y. The shift is caused by the dipolar field 
created on Rb sites by the field-induced magnetic mo- 




T(K) 

FIG. 8: (Color online)NMR spectra of RbFe(Mo0 4 ) 2 
recorded for a frequency of 90.6 MHz at various temperatures, 
H 1 C3 (top panel). Bottom panel shows the temperature 
dependence of the resonance fields. 

ments of the Fe 3+ ions. 

In the magnetically ordered phases, the dipolar fields 
of the Fe ions make the Rb ion positions non-equivalent, 
which results in splitting of all the spectral components 
into pairs of lines with distinct intensity. The more in- 
tense line shifts towards lower field, while the less in- 
tense line shifts towards higher field. The ratio of the 
intensities for these lines is approximately 2 to 1. The 
observed splitting can be used for a precise determina- 
tion of the Nccl temperature; it is also useful in choos- 
ing an appropriate model for the magnetic structure of 
RbFe(Mo04)2^- The values of the ordering tempera- 
tures, as inferred from the NMR results, are also collected 
in Fig. 1111 together with the magnetisation and specific 
heat points. 

In order to determine the temperature dependence of 
the order parameter, detailed NMR measurements were 
performed near the critical temperature for one of the 
spectral components. This component is seen at a fre- 
quency of 81 MHz in a field of 5.82 T, which is nearly 
in the centre of the collinear phase. The values of the 
resonance fields for a central line above and below its 
splitting point are plotted in Fig. [SJ 

Let us define here the order parameter in the collinear 
phase following Ref. 
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FIG. 9: Temperature dependence of the NMR resonance field 
B T es{T) and the magnetic moment M(T). fioH ~ 6 T, H 1 
C 3 . 



t ? =((M 1 +M 2 )/2)-M 3 )/2, 
where Mi, M2 and M3 are the magnetic moments of the 
three sublattices normalised per magnetic ion. The pa- 
rameter r\ is zero in the paramagnetic phase and nonzero 
in the ordered collinear phase. Let us also introduce a 
normalised full magnetic moment: 

/*= (Mi +M 2 +M 3 )/3. 

The temperature dependence of the order parameter 
r\ can be established by measuring the local field B e ff, 
which is created by the Fe 3+ ions on the Rb nuclei. Ex- 
perimentally, this field is defined as 

B ef f = B res - 1//7 - Av q /j, 
where Av q is the quadrupole shift of the resonance field. 

At any finite temperature the magnetic structure 
can be represented as a superposition of the ferromag- 
netic and collinear (UUD) states with the corresponding 
weights fi— 77/3 and 77, respectively. Therefore, the effec- 
tive field B e ff consists of two components, one of them 
being proportional to /1, and the other being proportional 
to rj(T). The order parameter is hence given by the ex- 
pression 

rj(T) ~ B eff (T) - B eff (T ) ■ M(T)/M(T ), 
where To is a fixed temperature above the ordering tem- 
perature Tff. From this expression, by using the indepen- 
dently measured temperature dependence of the mag- 
netic moment, one can determine the order parameter. 
Fig. [5] gives the temperature dependence of the NMR 
resonance field as well as of the magnetic moment M(T) 
measured in a field of ~ 6 T. Fig. 1101 shows the temper- 
ature dependence of the order parameter obtained us- 
ing the data of Fig. as described above. An absolute 
value of the order parameter is obtained by using the re- 
sults of Ref. 0, where the value of 77 was determined for 
T = 1.6 K and a magnetic field of 6 T. An abrupt change 
in the order parameter is clearly seen near the ordering 
temperature. 
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FIG. 10: Temperature dependence of the order parameter 77. 
The solid and open circles correspond to the values of 77 ob- 
tained from the resonance fields of 87 Rb residing in different 
positions with respect to the magnetic structure. The trian- 
gles correspond to the values of 77 in the paramagnetic phase. 
ftoH ~ 6 T, H 1 C 3 . 



C. Magnetic phase diagram 

Fig. II II summarises the positions of the phase bound- 
aries in RbFe(MoC>4)2 for a field applied in the basal 
plane of the crystal, determined by using different meth- 
ods. In weak fields, the transition temperature is well de- 
fined (with the accuracy better than 0.1 K) from the mag- 
netisation and specific heat measurements. The phase 
boundary between the paramagnetic and the collinear 
phase (P3 in Fig. fTT|) is also well defined by the peak 
in the specific heat and by the splitting in the NMR 
spectrum, while the magnetisation curve shows a smooth 
behaviour. The phase boundaries corresponding to the 
transition fields H C 2 and H c % are only visible in the mag- 
netisation curves, whereas the anomaly at H c \ is clearly 
seen in both the specific heat and the magnetisation 
curves. 

The dashed lines in Fig. ^Jrepresent the phase bound- 
aries obtained in our earlier paper—. Although the over- 
all agreement between the two sets of the results is good, 
a noticeable difference between the transition tempera- 
tures is most likely to be caused by the variations in 
sample quality, originating from difference in the sample 
preparation procedures. The specific heat peaks associ- 
ated with the phase transition from the paramagnetic to 
the collinear phase were much sharper for the new batch 
of samples, for which the NMR signal splitting also occurs 
over a narrower temperature interval. During the mea- 
surements with different samples from the same growth 
batch it has been noticed that by rigidly gluing the NMR 
sample to a sample holder it was possible to broaden sig- 
nificantly the phase transition region. In addition, the 
specific heat measurements (performed as usual for the 
heat-pulse method with the crystals attached rigidly to 
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FIG. 11: (Color online)Magnetic phase diagram of 
RbFe(Mo04)2 for H _L C3 determined by using different ex- 
perimental methods. Spin structures of magnetic layers are 
schematically presented by solid arrows and marked as A, 
B, C. Two-layer fragments of supposed 3D spin structures 
are shown by solid and dashed arrows for different magnetic 
phases marked as P1-P5, see text. 



a platform using a low-temperature grease) show that 
the peaks in C (T) measured for the PM to P3 transition 
in the thinner and therefore more stressed samples were 
about 0.2 K broader compared to the thicker and less 
stressed samples. It should be added, however, that the 
shape and the width of the peaks in C(T) curves were 
not sensitive to the stress for the other magnetic phases 
(above and below the phase with a collinear structure). 
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FIG. 12: Temperature dependence of the specific heat in dif- 
ferent fields (top panel) and spectra of the central quadrupole- 
split NMR line measured at 35 MHz for different temperatures 
(bottom panel). The field is applied along the C 3 axis of the 
RbFe(Mo04)2 single crystal. 



IV. MAGNETIC PHASE DIAGRAM FOR H || C* 3 

Fig. 1121 (upper panel) shows the temperature depen- 
dence of the specific heat measured in RbFe(MoC"4)2 for 
H || C 3 . The lower panel also shows the records of the 
central quadrupolar-split NMR line measured at 35 MHz 
for the same direction of the field. When the field is ap- 
plied along the C 3 -axis ( "hard" axis in terms of the mag- 
netic anisotropy), a single phase transition into a mag- 
netically ordered phase is observed via a A-like anomaly 
in C(T). This transition is seen as a sharp maximum 
in the C(H) curves (not shown) and also as a splitting 
of the NMR line. For H || C 3 , this splitting is nearly 
an order of magnitude smaller compared to the H _L C 3 
geometry. 

The H — T phase diagram for a field applied along the 
C 3 -axis is shown in Fig. The circles correspond to 
the phase boundary positions obtained from the specific 
heat measurements, while the square marks the phase 
transition observed by the NMR splitting. 



V. DISCUSSION 
A. Initial aspects 

Here we discuss the observed phase transitions in 
RbFe(Mo04)2 and its critical properties. Apart from 
the data presented in this paper, the discussion will be 
based on the available data on the microscopic structure 
of RbFe(Mo04)2 obtained from neutron diffraction^ and 
NMR 16 experiments, as well as on the established hier- 
archy of the magnetic interactions^: the intralayer ex- 
change, the easy-plane anisotropy and the interlayer ex- 
change. Besides these experimental facts we shall use 
the results of the calculation of the dipole energy, which 
can help to select between the proposed states with equal 
nearest neighbors exchange energy. 

It has to be noted that in zero field and at finite tem- 
perature for the case of an ideal two-dimensional system, 
long range order of the magnetic components is prohib- 
ited by the Mermin- Wagner theorem. For a triangu- 
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FIG. 13: The magnetic phase diagram of the RbFe(Mo04)2 
single crystal for H || C 3 . The dashed line corresponds to 
the relation /i H = 20.2 T x(l — T/3.9 K) ' 5 . The value 
of the experimental uncertainty in the position of the phase 
boundary is represented by the markers size. 



lar lattice with antiferromagnetic exchange interactions, 
however, apart from the Berezinskii-Kosterlitz-Thouless 
transition there should be an additional transition to a 
state with long-range order with the chirality parameter 
for the neighbouring triangles, the so-called "staggered 
helicity (vorticity)" stated. Monte Carlo simulations 
for this stated show that this phase is close to the state 
with long-range order of the nonzero average magnetic 
moments, as the spin-spin correlations decay in a power- 
law manner. Moreover, the characteristic relaxation time 
for the 120-degree three-sublattice structure could be rel- 
atively long and may exceed the dynamical range for var- 
ious experimental techniques. 

Computer simulations show the presence of chiral or- 
dering both for the 2D XY-model^ and 2D Heisenberg 
systemAii For the 2D XY-model, in addition to the 
paramagnetic phase, the H — T phase diagram contains 
the three ordered phases^ which have previously been 
identified in, for example, Ref. |8j These structures are 
depicted and labelled as A, B and C in Fig. ^J] 

In the presence of an external magnetic field, the sta- 
bilisation of long range order becomes inevitable^. Com- 
puter simulations for the Heisenberg spins 9 *i£ reveal qual- 
itatively similar results to the XY-model phase diagram, 
however, the exact positions of the phase boundaries are 
determined there far less accurately. Both models predict 
a characteristic non-monotonic variation of the ordering 
temperature with the magnetic field. The simulation re- 
sults for the XY-model are shown in Fig. ITU 

Weak interlayer exchange coupling should, in the gen- 
eral case, facilitate the appearance of 3D long range or- 
der. Given the results for the 2D models discussed above 
and the proximity of the ground states for these systems, 
one can speculate that the temperature of the 3D order- 
ing should be close to the critical temperature of the 2D 



transition. In order to corroborate this conjecture, ad- 
ditional Monte Carlo simulations for the system with a 
finite interlayer exchange and magnetic anisotropy would 
be required. 

In order to analyse the real structures in 3D crys- 
tals one has to take into account the interlayer ex- 
change coupling and to consider various possibilities for 
the alignment of the spins in the neighbouring planes. 
Theoretical analysis of the different allowed structures 
for an antiferromagnetic interlayer exchange has been 
performed in Ref. 6. The proposed magnetic phases 
have then been used as the models for RbFe(MoC>4)2 in 
Ref. E3 Fig. UU shows the spin configurations of the 
ordered phases, which have been selected from the pro- 
posed \\s& on the basis of matching the NMR-detected 
commensurate-incommensurate transition 16 , magnetisa- 
tion plateau, and possible periodicity along C 3 direction, 
which is discussed below. Solid and dashed arrows in 
Fig. 1111 represent the orientations of the magnetic mo- 
ments of the Fe 3+ ions from neighboring planes for these 
selected configurations. 

Let us consider the procedure of the selection of the 
proposed structures. One should note that for the spin 
structures with two parallel sublattice magnetisations 
within a layer (e.g. structures suggested for phases 
P2,P3,P4), an interlayer exchange does not define the 
magnetic periodicity along the C 3 -axis, which could 
therefore be any integer multiple of the lattice periods. 
This point is illustrated in Fig. 1151 where the two struc- 
tures with the periodicity along the C 3 -axis equal 2c and 
3c, contain the same number of various angle combina- 
tions for the neighbouring spins. These two structures 
thus possess the same exchange energy. This periodicity 
must be due to some other interactions, that are weaker 
when compared to the exchange coupling. Among the 
possible candidates for the interactions which select a pe- 
riod along C 3 axis are the dipole-dipole interactions and 
further neighbour exchange interactions along the C - 
axis. The direct calculation of the dipole energy shows 
that it is lower for the period of 3c. This gain in the 
dipole energy equals approximately 0.4 mK per Fe-ion 
for UUD structure. Therefore we suppose that the struc- 
tures suggested for phases P2, P3, P4 have the period 3c. 
The next-nearest antiferromagnetic exchange interaction 
in the C 3 direction would also result in the preference of 
period 3c. 

Fig. II II shows results of the various experimental tech- 
niques that reveal the presence of five ordered phases. 
The phase transition, corresponding to the field H C 3, 
which is marked in the 4ff curves in Fig. [21 is an 
incommensurate-commensurate transition indicated by 
the characteristic change of the NMR spectrum^. We 
suppose that the incommensurate structure below H C 3 
corresponds, as usual, to a nearly antiparallel orienta- 
tion of spins from neighboring layers (structure PI), this 
structure being rotated around the C 3 axis by a small 
angle from layer to layer. In the range above H c3 the 
commensurate structure takes place for the phase P2. In 



9 




XY-model 



Experiment a ^a^«* 

I^a-a^ • ^*«*-aa a .^| 



0.0 



0.2 



0.4 



0.6 



T/JS 



FIG. 14: (Color online) Experimental (open symbols) and 
theoretical- (solid symbols) temperature dependence of the 
specific heat in RbFe(Mo04)2 plotted against the normalised 
temperature in zero field. The inset shows the H — T phase 
diagram for H _L C3, where the open symbols represent the 
experimental results as in Fig. 1111 solid symbols are the results 
of computer simulations- with J = 1.23 K. The solid lines in 
the inset are drawn as a guide to the eye. 



Ref. 0, two possible variants were proposed for the spin 
structures in phases PI and P2 based on the results of 
theoretical considerations 6 ^. We selected for the phase 
P2, the structure which allows the period of 3c, favorable 
for dipole-dipole energy. For the phase PI the structure, 
compatible with period 2c is suggested, since for the typ- 
ical incommensurate structures in antiferromagnets the 
orientation of spins in neighboring layers should be close 
to antiparallel one, like at the period value 2c. 

In the phase P3, i.e. in the plateau range, the collinear 
structure within a layer and the periodic structure along 

-axis, probably with the period 3c, is realised. The 
structures of P4 or P5 type may be proposed for the 
high-field range between H C 2 and saturation. We have 
here no NMR data which could distinguish between com- 
mensurate and incommensurate states, thus we can not 
determine the nature of the P4-P5 boundary, detected 
by small dM/dH singularity in Rcf 12 «. Nevertheless, if 
we extrapolate the period value of 3c from P3 to P4 , 
the structure which is compatible with period 3c should 
be selected, as shown in Fig. In this case the transi- 
tion to the second high-field spin structure of the phase 
P5 should change the period to 2c or again result in the 
incommensurate modulation. 

The collinear structure of the P3 phase may be taken 
as determined most precisely, as it corresponds to the 
maximum value of the NMR shift 16 , which is maintained 
in the entire field range H c i < H < H C 2. The spin 
structures of other phases are suggested as probable con- 
figurations. 



B. 2D-3D crossover of the specific heat 
temperature dependence 

The overall shape of the phase diagram, as well as the 
characteristic widening of the collinear phase range with 
increasing temperature predicted in Ref. H points to the 
two-dimensional nature of magnetism in RbFe(Mo04)2- 
Here, we show that the observed change-over in the tem- 
perature dependence of the specific heat from the low- 
temperature T 3 law to the intermediate temperature T 2 
law (as described in section ITlI Aj) is also indicative of the 
two-dimensional disposition of this magnetic system. 

The temperature dependence of an antiferromagnet for 
T < 0.6 Tn is usually determined by thermally activated 
magnons. The spectrum of acoustic magnons in antifer- 
romagnetic RbFe(Mo04)2 consists of one gapless mode 
and two modes with a gap of 90 GHz (4.5 K)^ At such 
low temperatures, the variation of the gap's magnitude 
should be small, as the order parameter for T < 0.6 Tn 
is reduced by not more than 30%. The main contribu- 
tion to the specific heat in this temperature region comes 
from the magnons of the gapless branch. At the lowest 
temperature the wave-numbers of the thermal magnons 
are small and lie far away from the Brillouin zone bound- 
ary. Therefore one can expect for C(T) a typical three- 
dimensional behaviour, that is C(T) ~ T 3 . 

Due to a large difference between the inter- and the 
intra-layer exchange interactions, the magnon dispersion 
relation is highly anisotropic - the magnons propagat- 
ing along the C 3 -axis are less dispersive than those with 
a wave-vector confined to the triangular plane. Con- 
sequently, at a certain temperature, the energy of the 
thermal magnons travelling along the C 3 -axis reaches its 
maximum value, while the size of the A;-space region cor- 
responding to magnons with k _L C 3 will continue to 
expand. For these temperatures one can expect a typi- 
cal two-dimensional behaviour, that is C(T) ~ T 2 . In- 
evitably, the change-over from a cubic to a quadratic law 
should occur at the temperature dependent on interlayer 
exchange energy, or, more precisely, on the maximum en- 
ergy of the fc||C 3 magnons of the lowest gapless branch. 
This is the energy of the lowest "exchange" -mode of the 
antifcrromagnetic resonance spectrum. According to Ref. 
IT2I the gap of the lowest exchange mode is about 30 GHz 
(1.5 K), which agrees well with the observed crossover 
temperature. Thus the 2D-3D crossover in the temper- 
ature dependence of specific heat confirms the domina- 
tion of two-dimensional features in the thermodynamic 
behaviour of RbFe(Mo04)2- 



C. Critical properties of the magnetic ordering 
transition for H _L C 3 

Given the two-dimensional nature of the magnetic 
properties of RbFe(Mo04)2, let us compare the ob- 
served critical behaviour with the results of the classi- 
cal Monte Carlo simulations^ performed for the 2D XY- 
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n+2 

FIG. 15: Schematic representation of the magnetic structures 
with different periodicity along the C 3 -axis: 2c on the left (a) 
and 3c on the right (b). Solid and dashed arrows correspond 
to the direction of the magnetic moments in the neighbouring 
planes. Numbers 1, 1' and 1" mark the magnetic moments of 
the Fe 3+ ions located above one another in the neighbouring 
n th , n + 1 st and n + 2 nd magnetic layers. 



model. The Neel temperature in this simulation is de- 
fined as 0.5 JS 2 . Using the value of J = 1.2±0.1 K deter- 
mined from the susceptibility and saturation magnetisa- 
tion measurements 1 ^*!—, the predicted ordering tempera- 
ture is 3.75±0.3 K, which is in a good agreement with the 
value of Tn — 3.85 ± 0.05 K determined experimentally. 
The ordering temperature for the 2D Heisenberg antifer- 
romagnet on a triangular lattice, on the other hand, is 
defined as 0.35 JS 2 ^ which unexpectedly gives less sat- 
isfactory agreement with the experimental T^. 

Fig. 1141 shows the comparison between the tempera- 
ture dependence of the specific heat measured experimen- 
tally and obtained from the Monte Carlo simulations*. It 
also shows a comparison of the phase boundary positions 
from the experimental H — T diagram (see Fig. Illfl and 
the simulations 4 for H _L C 3 . For the theoretical curve, 
the value of the exchange parameter quoted above has 
been used. The overall agreement is satisfactory, as the 
two phase diagrams have very similar shapes, apart from 
the low- field region. For low fields, however, the dis- 
tinction is fundamental: the theoretical model predicts a 
phase transition from paramagnetic to chiral structures 
only through a collinear phase, while experimentally a 
single phase transition is observed for all fields below 2 T. 

The measured C(T) and 77 (T) curves allow one to de- 
termine the specific heat and the order parameter critical 
indices for the phase transition into the collinear phase. 
Corresponding experimental data are presented in Figs. 
1101 and 1161 The experimentally determined indices are 
a = 0.40 ± 0.03 for the specific heat and j3 = 0.13 ± 0.02 



10 




T(K) 

FIG. 16: Temperature dependence of the specific heat for a 
field of (ioH — 5 T applied in the easy-plane of RbFe(Mo04)2. 
The inset shows C vs 1 — T/T c \ curve in the vicinity of the 
phase transition on a log- log scale. 

for the order parameter. They agree well with a = 1/3 
and (3 = 1/9 obtained from the simulations*. 

It can be concluded then that the ordering tempera- 
ture, the magnetic phase diagram and the critical indices 
of the specific heat and the order parameter demonstrate 
a satisfactory quantitative agreement with the 2D XY- 
model predictions for zero field and also for regions of 
the field (H ± C 3 ) where the collinear phase is stable. 
Most likely this agreement can be explained by the large 
value of the energy gap for the magnon branch, which 
corresponds to a deviation of the spins out of the basal 
plane. 

This gap has been estimated as 90 GHz (or 4.5 K)*£. 
Therefore, for all the magnetically ordered phases, the 
thermodynamic properties are defined by gapless and 
low-energy magnons, corresponding to spin fluctuations 
within the easy plane. This, in turn, justifies the appli- 
cability of the XY model. 



D. Phase diagram for H || C 3 

The easy-plane type of magnetic anisotropy in 
RbFe(Mo04)2 is confirmed by the observation of the 
magnetisation plateau around M s /3 for the magnetic 
field applied in the ab plane and by the linear rise of 
the M(H) curves for the field parallel to the C 3 -axis. 
For an easy-plane anisotropy, the magnetic moment of 
the system is generated for H || C 3 by the even tilts of 
all the subblattices in an umbrella-like structure towards 
the field (see Fig. Il"3"|) . For such a system, the application 
of an external magnetic field should result in a decrease 
of the ordering temperature. Moreover, in the vicinity of 
Tn{H = 0), the reduction should be proportional to the 
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square of the field, similarly to an ordinary two-sublattice 
antiferromagnefeiS. The dashed line in Fig. ^|shows such 
a dependence, H c = o(l - T/T c ) 5 , with a = 20.2 T and 
T c = 3.9 K. For the case of the uniformly tilted spins 
described here, splitting in the NMR spectra should be 
absent, but it is observed. The origin of this small but 
clearly visible splitting remains unknown at the present 
time. 



VI. CONCLUSIONS 

To summarise, we have performed a systematic and 
comprehensive study of the thermodynamic and mag- 
netic properties of the quasi two-dimensional TLAM 
RbFe(Mo04)2, as well as measurements of the magnetic 
order parameter for this compound. A crossover from a 
T 2 law to a T 3 law is observed for the temperature depen- 



dence of the specific heat at low temperatures. The phase 
diagram and critical properties demonstrate a surpris- 
ingly good quantitative agreement with the predictions 
for the XY-model despite the fact that the easy-plane 
anisotropy in RbFe(Mo04)2 does not dominate the main 
exchange interaction. 
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